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We report the results of magnetization and specific heat measurements on Ba3NiSb209, which 
is a quasi- two-dimensional spin-1 triangular- lattice antiferromagnet. We observed a nonclassical 
magnetization plateau at one-third of the saturation magnetization that is driven by spin 
frustration and quantum fluctuation. Exact diagonalization for a 21-site rhombic cluster was 
performed to analyze the magnetization process. Experimental and calculated results agree 
well. 

KEYWORDS: Ba3NiSb209, triangular-lattice antiferromagnet, magnetization plateau 



I 

o 
o 



> 

o 

in 
o 



Triangular-lattice antiferromagnet (TLAF), which is 
one of the simplest frustrated systems, exhibits a rich 
variety of physics. ^"^^ A symbolic quantum effect pre- 
dicted for small-spin Heisenberg TLAF is that magneti- 
zation is quantized at one-third of the saturation mag- 
netization Ms in magnetic fields For classical spin, 
the equilibrium condition in a magnetic field H is given 
by Si -\- S2 + S3 = g^^H I (3 J) with the sublattice spins 
Si . Because the number of the equations is smaller than 
the number of parameters to determine the spin configu- 
ration, the stable state is not uniquely determined; thus, 
the ground state is infinitely degenerate. This classical 
degeneracy can be lifted by the quantum fluctuation that 
is remarkable for small spin, and a specific spin state is 
selected as the ground state. Consequently, the up-up- 
down state is stabilized in a finite field range, which can 
be observed as a magnetization plateau at \Ms 

In experiments, the quantum magnetization plateau 
at \Ms was actually observed in Cs2CuBr4 ^n^i^) -^^j^gj-g 
Cu^ ions having spin-i form a spatially anisotropic tri- 
angular lattice with J'/J = 0.74. In Cs2CuBr4, an addi- 
tional quantum plateau was also observed at |Ms 
Furthermore, precise thermodynamic measurements re- 
vealed the cascade of field-induced quantum phase tran- 
sitions.^^-' These unusual observations should be at- 
tributed to the spatially anisotropic triangular lattice 
and the Dzyaloshinsky-Moriya interaction but the 

overall explanation is still an open question. 

The motivation of this study comes from a question on 
whether the ground states for Heisenberg TLAFs with 
integer and half-integer spins are qualitatively the same. 
It is known that in the case of kagome-lattice Heisen- 
berg antiferromagnet that is closely related to Heisen- 
berg TLAF, the ground states for spin-i and spin-1 are 
qualitatively different For spin-1 Heisenberg TLAF, 
the ground states in zero and finite magnetic fields have 
not been well understood. NiGa2S4 is known as spin- 
1 Heisenberg TLAF with incommensurate spin corre- 
lation .24-26) jj-^ -(-j^jg nickel compound, no clear magne- 




Fig. 1. (Color) Crystal structure of Ba3NiSb209. Dotted lines 
denote the chemical unit cell. 



tization plateau has been observed, although the mag- 
netization shows a tiny undulate anomaly at around 
\Ms To investigate the magnetization process of 
spin-1 Heisenberg TLAF, we have performed high-field 
magnetization measurements on Ba3NiSb2 0g. As shown 
below, a nonclassical magnetization plateau was observed 
at \Ms. 

Ba3NiSb2 09 crystallizes in a high symmetric hexago- 
nal structure, P&^/mmc^'^^''^'^^ which is the same as the 
hexagonal BaTiOs structure (see Fig. 1). The structure is 
composed of single NiOg octahedra and face-connected 
Sb2 09 double octahedra, which are linked sharing cor- 
ners. Magnetic Ni^+ ions form layers of triangular lat- 
tice parallel to the c plane, which are separated by the 
nonmagnetic layer of Sb209 double octahedra. There- 
fore, the interlayer exchange interaction is expected to 
be much smaller than the intralayer exchange interac- 
tion. Ba3NiSb209 is considered to be one of the best 
candidates of spin-1 Heisenberg TLAF. However, owing 
to the weak interlayer exchange interaction, Ba3NiSb209 
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Fig. 2. (Color) Specific heat of Ba3NiSb209 measured at zero 
magnetic field. Inset shows an expansion at around 13 K. 



undergoes antiferromagnetic ordering at Tn — 13.5 K ."^"^ 
In this study, we also measured the specific heat to in- 
vestigate the nature of the magnetic ordering. As shown 
below, successive phase transitions were observed. 

Ba3NiSb20g powder was prepared via a chemical re- 
action SBaCOg -I- NiO + SbsOs — > BagNiSbaOg -f- 
3C02- Reagent-grade materials were mixed in stoichio- 
metric quantities, and calcined at 1200 °C for 30 h in 
air. Ba3NiSb209 was sintered at 1200 1600 °C for more 
than 30 h after being pressed into a pellet. The sam- 
ple obtained was examined by X-ray powder diffraction. 
A tiny peak due to nonmagnetic impurity BaSb206 of 
~ 2 wt% was observed. 

The specific heat of Ba3NiSb209 was measured down 
to 1.8 K using a physical property measurement system 
(Quantum Design PPMS) by the relaxation method. 
High-field magnetization measurement of up to 53 T 
was performed using an induction method with a mul- 
tilayer pulse magnet at the Institute for Solid State 
Physics, University of Tokyo. We also measured the tem- 
perature dependence of magnetic susceptibility and con- 
firmed that the result was the same as that reported by 
Doi et al. .^o) 

Figure 2 shows low-temperature specific heat mea- 
sured at zero-magnetic field. As shown in the inset of 
Fig. 2, Ba3NiSb209 undergoes two magnetic phase tran- 
sitions at Tni = 13.5 K and Tn2 — 13.0 K. In Heisenberg- 
like TLAF, the successive phase transition arises when 
the magnetic anisotropy is of easy-axis type, while a 
single transition occurs for the easy-plane anisotropy. 
The transition scenario for the easy-axis anisotropy is as 
follows: with decreasing temperature, the c axis compo- 
nent of spins orders first at Tni and then the ab com- 
ponents order at Tjq2- Consequently, below Tn2, spins 
form a triangular structure in a plane parallel to the c 
axis. The reduced temperature range of the intermedi- 
ate phase (Tni — Tn2)/?ni is determined from the ra- 
tio of the easy-axis anisotropy to the exchange interac- 




Fig. 3. (Color) Magnetization curve, dM/dH vs H and dM/dH 
vs M of Ba3NiSb209 measured at 1.3 K. The data for dM/dH vs 
M is shifted upward by 1 X 10~2/iB/(Ni^+ T). A horizontal bar 
denotes the field range of the magnetization plateau calculated 
by exact diagonalization. 



tion.31.32) Tj^g 

successive phase transitions with a very 
narrow intermediate phase in Ba3NiSb2 09 indicate that 
the magnetic anisotropy is of the easy-axis type and is 
much smaller than the exchange interaction. The origin 
of the anisotropy should be a single-ion anisotropy of the 
form T>(5^)2. 

Figure 3 shows the magnetization curve of 
Ba3NiSb209 powder measured at 1.3 K and the 
derivative susceptibility dM/dH vs magnetic field H. A 
magnetization plateau is observed at M ~ 0.8 /iB/Ni^+. 
Because the g factor of Ba3NiSb2 09 powder is 
g~2.3,'^^-' the magnetization at the plateau corresponds 
to one-third of the saturation magnetization Mg. The 
magnetization anomalies at the lower and higher edge 
fields are smeared. In general, the g factors for H\\c and 
H Lc are slightly different and the field ranges of the 
plateau for these field directions are also different owing 
to the weak anisotropy energy. These facts lead to the 
smearing of the edge fields of the plateau in a powder 
sample. 

In the classical Heisenberg TLAF, stable spin state 
is not uniquely determined in a magnetic field due to 
spin frustration. However, small easy-axis anisotropy 
and thermal and quantum fluctuations can stabilize a 
coUinear up-up-down state in a finite field range. In the 
case of classical Heisenberg TLAF with the easy-axis 
anisotropy along the c axis, the ^-plateau appears for 
H \\c, while not for H _L c.^^^^^^ Thus, for a powder sam- 
ple of such classical system, dM/dH at the plateau range 
should be about two-third of dM/dH for H smaller than 
the lower edge field. In the present system, the mini- 
mum of dM/dH at the plateau range is smaller than 
half of dM/dH for <30T. This implies that the i- 
plateau appears irrespective of external field direction. 
In classical Heisenberg TLAF, thermal fluctuation stabi- 
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Fig. 4. (Color) (a) Magnetization curves and (b) derivative sus- 
ceptibility vs normalized magnetization M/Ms obtained by exact 
diagonalizations for 9-, 12-, and 21-site rhombic clusters. 



lizes the up-up-down state at finite temperatures, so tliat 
a plateaulike magnetization anomaly occurs in the mag- 
netization curve. "^^^ The magnetic field range of the ther- 
mally stabilized collinear state decreases with decreasing 
temperature and becomes zero at T = 0. The i-plateau 
in Ba3NiSb209 has a magnetic field range of about 10 T 
even at 1.3 K (<;Tn2). Thus, we can deduce that the 
i-plateau in Ba3NiSb209 arises from the quantum fluc- 
tuation .^"^^ However, the magnetization at the plateau 
range is not completely flat. This should be ascribed to 
the increase in the sample temperature due to the mag- 
netocaloric effect and the susceptibility of impurities. 

Following the method described by Troyer et 
al.^^^ for one-dimensional gapped excitations, we 
can derive the low-temperature molar susceptibility 
X = iVAnmg^^B 6xp(— /3A)/(27ra), for two-dimensional 
gapped excitations with dispersion e{k)—A + ak^. 
Here, is the number of gapped modes and 

rim = 2 in the present case. When the external field 
is one-third of the saturation field Hs, the coeffi- 
cient a is given by a = (7/8)J.^^ If the susceptibility 
dM/diJ=8xlO-3^B/(Ni^"^T) at the center field of the 
^-plateau arises only from the finite temperature effect, 
the spin temperature is estimated to be r~2.1K using 
A/kB ^ 4.5 K39) and J/ks = 19 K, which is given below. 

Figure 4(a) shows magnetization curve of spin-1 
Heisenberg TLAF calculated by exact diagonalization for 



9-, 12- and 21-site rhombic clusters. The magnetiza- 
tion curve for spin-1 Heisenberg TLAF using exact di- 
agonalization has already been reported by Lauchli et 
al..^"^^ However, their calculations for iV> 21 are fimited 
to a higher-field part of the magnetization. Note that our 
result for TV = 21 covers all the magnetization values. In 
Fig. 4(b), we also present derivative susceptibility x de- 
fined as 

_i _ E{N, M + 1) - 2E{N, M) -f E{N, M - 1) 
^ " l/M, ' ^' 

where E{N^ M) is the lowest energy for a given z com- 
ponent of the total spin S^^^^ = M for a given cluster size 
N . It is easier to detect the critical behavior around 
the plateau in x rather than only in the magnetization 
curve. Our results in Figs. 4(a) and (b) reveal a mag- 
netization plateau at \Ms, although the magnetic field 
range of the plateau is fairly small as compared with the 
spin-i case.'''^''*^^ 

As shown in Fig. 4, the calculated susceptibility in- 
creases with increasing magnetization and shows a max- 
imum just below M = ^Ms- Just at M = |-Ms, x discon- 
tinuously becomes small. This behavior is also observed 
in dM/dH vs M shown in Fig. 3. For M > ^Ms, the 
calculated magnetization and susceptibility exhibit no 
anomaly suggestive of quantum phase transition up to 
saturation. 

The lower and higher edge fields tend to converge to 
Hci =2. 85 J and i?c2 = 3. 72 J, respectively, with increas- 
ing cluster size. A horizontal bar in Fig. 3 denotes the 
calculated field range of the i-plateau. Here, we chose 
the center field of the ^-plateau as {Hd + Hc2)/2 = 41 T. 
The experimental field range of the plateau is roughly es- 
timated as ^ 60% of the calculated field range. From the 
center field of the plateau of 41 T and 5^2.3, the ex- 
change constant is evaluated as J/fce — 19 K. This value 
of J is approximately | times as large as J/ fee — 32 K 
that is obtained from the Weiss constant 8 = — 128 K 
and the molecular field theory. 

For the spin-1 case, the field range of the ^- 
magnetization-plateau normalized by the saturation field 
is calculated as {Hc2 — Hci)/Hs = 0.097, while for the 
spin-i case, the normalized plateau range is O.lTl.^"'*'^'' 
The ^-magnetization-plateau for the spin-1 case is con- 
siderably suppressed as compared with that for the the 
spin-i case. The field range of the i-plateau and the 
dM/dH vs M observed in Ba3NiSb2 09 are consistent 
with those calculated for the spin-1 Heisenberg TLAF. 
This indicates that the magnetization plateau observed 
in Ba3NiSb2 09 arises from the quantum effect charac- 
teristic of the spin-1 Heisenberg TLAF. 

For TLAF with S >1, the magnetization plateau at 
iMg can be stabilized with the help of the biquadratic 
exchange interaction of the form —Kij{Si-Sj)'^ with 
/Cy > 0,42,44) ^Yiich mainly arises from the spin-lattice 
coupHng as discussed by Penc et aZ..*^^ The biquadratic 
exchange interaction may be responsible for the magneti- 
zation plateau observed for H _L c in RbFe(Mo04)2,^^"''^' 
in which the normalized plateau range is as large as 0.10 
despite spin-|. In Ba3NiSb2 09, the contribution of the 
biquadratic exchange interaction should be minimal, be- 
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cause the field range of tfie magnetization plateau can 
be explained quantitatively within the bilinear exchange 
interaction. 

As shown in Fig. 4(b), the calculated susceptibility ex- 
hibits a sharp dip at M = 0, almost irrespective of clus- 
ter size N. No such anomaly is observed for the spin- 
I triangular- and kagome-lattice antiferromagnets with- 
out a spin gap.'*^''**^ At present, the origin of the sus- 
ceptibility anomaly at M = is unclear. On the other 
hand, the susceptibility observed at small magnetization 
in Ba3NiSb209 is finite because of magnetic ordering due 
to the weak interlayer exchange interaction. 

In conclusion, we have presented the results of spe- 
cific heat and high-field magnetization measurements on 
Ba3NiSb2 09 that is described as a spin-1 Heisenberg 
TLAF. The present system exhibits two phase transi- 
tions Tni = 13.5K and Tn2 = 13.0K due to weak easy- 
axis anisotropy. We observed a magnetization plateau 
at at temperatures below rx2. We performed ex- 

act diagonalization for rhombic spin chistcrs with up 
to 21-sitcs to analyze the magnetization process. The 
calculated results are in agreement with experimental 
observations. The magnetization plateau observed in 
Ba3NiSb2 09 arises from the interplay between spin frus- 
tration and the quantum effect characteristic of small- 
spin Heisenberg TLAF. 
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